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ABSTRACT 

The planet GJ 1214b is the second known super-Earth with a measured mass and radius. Orbiting 
a quiet M-star, it receives considerably less mass-loss driving X-ray and UV radiation than CoRoT- 
7b, so that the interior may be quite dissimilar in composition, including the possibility of a large 
fraction of water. We model the interior of GJ 1214b assuming a two-layer (envelope+rock core) 
structure where the envelope material is either H/He, pure water, or a mixture of H/He and H2O. 
Within this framework we perform models of the thermal evolution and contraction of the planet. 
We discuss possible compositions that are consistent with M p = 6.55 M®, R p = 2.678 i?©, an age 
r = 3 — 10 Gyr, and the irradiation level of the atmosphere. These conditions require that if water 
exists in the interior, it must remain in a fluid state, with important consequences for magnetic field 
generation. These conditions also require the atmosphere to have a deep isothermal region extending 
down to 80—800 bar, depending on composition. Our results bolster the suggestion of a metal-enriched 
H/He atmosphere for the planet, as we find water-world models that lack an H/He atmosphere to 
require an implausibly large water-to-rock ratio of more than 6:1. We instead favor a H/He/H20 
envelope with high water mass fraction (~ 0.5—0.85), similar to recent models of the deep envelope 
of Uranus and Neptune. Even with these high water mass fractions in the H/Hc envelope, generally 
the bulk composition of the planet can have subsolar water:rock ratios. Dry, water-enriched, and pure 
water envelope models differ to an observationally significant level in their tidal Love numbers of 
respectively ~ 0.018, ~ 0.15, and ~ 0.7. 

Subject headings: planets and satellites: general — planets and satellites: individual(GJ 1214b) 



1. INTRODUCTION 

Among the nearly 500 detected planet candidates, and 
in particular among the <~ 80 planets for which both 
the mass and radius have been determined, our Earth 
is the only planet that is known to harbor liquid wa- 
ter oceans on a solid surface crust. These conditions 
have proven favorable for the development of life forms. 
With the discove ry of CoRoT-7b (M p = 4.8 ± 0.8 M e ; 
Leeer et alJl200l and G J 1214b (M p = 6.55 ± 0.98 M e ; 
Charbonneau et aLll2009[) , the search for other habitable 
worlds has just recently passed the milestone of discov- 
ering transiting cxtrasolar planets in the 2-10 M® super- 
Earth mass regime. 

Higher-mass planets such as Uranus (M p = 14.5 M®) 
are predicted by interior models to retain a H/Hc- 
rich atmosphere (jHubbard et al.l ["19951) whose size may 
vary depending on M p , R p and temperature. Pres- 
sures at the bottom of this envelope are high, and 
may reach 1 to 1000 kbar if the planet is of Neptune- 
size (INettelmann et aiT1l2010D or even several Mbar if of 
Saturn-size (|Guillotill999T before a presumably solid core 
is reached. Lower-mass objects on the other hand such 
as Mars or Ganymede have been observed not to retain a 
thick enough atmosphere that could prevent the planet's 
surface from cooling below the freezing point of water. 

CoRoT-7b -orbiting a Sun-like star at short orbital 
distance a p = 0.017 AU- falls into the intermediate 
super-Earth mass regime; the stellar extreme ultravio- 
let (EUV) flux it receives is so strong that its current 
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atmosphere is either a tiny remnant of an initially mas- 
sive gaseous envelope, or hot evaporating core material 
at an equilibrium temperature T K „ = 1 800 — 2600 K 
(jValencia et al.l [20101 1 Jackson et all [20Toh . In contrast, 
GJ 1214b (a p = 0.01 44AU) orbits an M star of ~ 3x 10~ 3 
smaller luminosity ([Charbonneau et al.l 120091) translat- 
ing into a planet-average T cq < 555 K, only a factor of 
two higher than that of Earth. Hence the discovery of 
GJ 1214b manifests an important step toward a detection 
o f an extrasolar ocean pla net . 

iRogers fc Seagerl (|2010f) investigated the response of 
interior models to the uncertainties in M p , and in- 
trinsic temperature Ti, lt and showed that GJ 1214b might 
have a (10~ 4 — 0.068) x M p thick H/He atmosphere, or 
else an outgasscd H atmosphere, or a water envelope atop 
a silicate-iron core with an ice:rock (I:R) ratio of 0.06 to 

CO. 

In this paper we adopt the fiducial {M p , R p } values 
and investigate how the unknown temperatures of the 
deep interior can be constrained by thermal evolution 
calculations. Our models are two-layer models with one 
homogeneous envelope overlying a rock core. We take 
into account mass loss during evolution and explore how 
that affects the possible mass of an outer H/He layer 
(ij l3.ll) . In § 13.21 we consider pure water atmospheres and 
ask whether condensation or even solidification of wa- 
ter could then have occured within 10 billion years of 
cooling. In line with recent transmission spectrum mea- 
surements that i ndicate 70% or m ore water by mass in 
the atmosphere ([Bean etaDEOlS), we vary in § 13.31 the 
envelope water mass fraction between 50 and 100% and 
suggest plausible models with about 0.2 x solar I:R ratio. 
In §2] we discuss our model assumptions and propose to 
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discriminate between our three classes of models (dry, 
water, water-rich envelope) by observationally determin- 
ing the Love number fc 2 and the mean molecular weight 
of the atmosphere. Our method of modeling this planet 
is explained in §[2j where we describe the irradiated at- 
mosphere grid (§[2J}, applied equations of state (§ 12. 2j) . 
mass loss ( §I2.3|) . and the calculation of structure and 
evolution (§' |2~4[) . 

2. METHODS 

In this section we describe the four components that 
our interior models rely on: the model atmosphere grid, 
the equations of state used, structure assumptions, and 
the thermal evolution to the present state. The evolu- 
tionary models have some similarities to what has pre- 
viou sly been applied to hot Neptunes and h ot Jupitcrs 
(e.g. iFortnev et all 120071 : iBaraffe et all 120081 ). but with 
additional complications due to radiogenic heating. As 
the planet's i nterior cools, the ex ternal radiative zone 
grows deeper (|Guillot et al.|[2006l) . reaching a depth of 
up to several hundreds of bars. The transition pressure 
of the atmosphere from radiative to adiabatic, at the 
current time, P a d(to), is a quantity we aim to constrain 
with our evolution model. We also perform explorations 
of the planet's structure as a function of Pad{to), to in- 
vestigate the full range of hotter, higher entropy inte- 
riors (lower P a d(to)) and cooler, lower entropy interiors 
(higher P a d(to)) that may be possible today. 

2.1. The model atmosphere 

Under the assumption the deep envelope layers con- 
vect efficiently, it is the radiative atmosphere atop the 
convective region that serves as the bottleneck for in- 
terior cooling, just as i n Neptune, Ur anus, and giant 
planets generally (e.g., IHubbard! 11977ft . Observations 
of GJ1214b's atmosphere are consistent with a water- 
dominated composition as w ell as with a H/H e atmo- 
sphere with clouds or hazes (jBean eEaDHQlO). Given 
our current ignorance of the composition of the atmo- 
sphere, we consider the two likely end-member cases, 
either a H/He-dominatcd atmosphere, or a pure steam 
atmosphere, and find tha t the opacity is dominated b y 
water vapor in either case ([Miller-Ricci fe Fortnevl[2010l ) . 
For planetary structure, we assume chemical equilibrium 
in the atmosphere, thereby ignoring possible alterations 
of the T-P profile through photoionization. 

For planetary evolution, a grid of model atmospheres 
is generally u sed as the upper bou ndary condition, see, 
for instance, IFortnev et al.l ([20071 ). These grids relate 
the specific entropy (s) of the convective interior, sur- 
face gravity (g) of the planet, and the intrinsic effec- 
tive temperature (Tint) from the interior. We have com- 
puted such a grid from T nlt = 175 K down to 30 K, 
with the correct limiting behavior down to K (an ex- 
hausted interior) across surface gravities from 100 to 1500 
cm s~ 2 . The grid is computed at 50 x solar metallic- 
ity, under the assumption of complete re distribution of 
absorbed stellar flux (meaning f=l /4, see IFortnev et all 
I2007t iMiller-Ricci fc Fortnevll2010( ). and no clouds. Sim- 
ilar models, which describe the te chnique in more detail, 
are found in IFortnev et al.l (120071 ) . The opacity database 
is described in iFreedman et al.l ([2008D . We note that 
this very high metallicity is realistic given that Nep- 



tune and Uranus are 30 -60 times solar in carbon (see 
IGuillot fc Gautien (120091) . for a review). 

We use this grid for all evolution calculations of our 
models, whether they posses thin H/Hc - atmospheres or 
pure steam atmospheres. This is certainly a broad brush 
treatment for a wide range of possible atmospheres, but 
given our current ignorance regarding the planet's atmo- 
sphere, we feel our treatment is justified. The importance 
and utility of the coupled model atmosphere/interior 
cooling calculation is that it allows us to estimate T lnt 
and Pad as a function of time. Fo r instance, in the recent 
work of IRogers" fc Seagerl pOlOf ). the value of T int was 
not calcul ated, but was extrap olated from evolutionary 
models of Baraffe ct al.l (|2008| ). for higher mass objects. 
Generally, we find a 15 K lower Tint, meaning a colder 
interior, than IRogers fc Seagerl (l2010t) used. While T int 
may change from 175 down to 30 K during evolution, we 
find that the effective temperature T e g remains nearly 
constant within 562 to 557 K, which is close to the zero- 
albedo, planet-average equilibrium temperature of 555 K. 

2.2. Equations of state 

Metal-rictH planets such as super-Earths are gener- 
ally suspected to harbour a variety of materials. We 
aim to represent this variety in a simplified manner by 
confining silicates and iron into a 'rocky' core, and H, 
He, and water to an envelope. F or core material we 
use th e P — p relation for rocks by IHubbard fc Marlevl 
(|1989t) which describes an adiabatic mixture around 
10 4 K of 38% Si0 2 , 25% MgO, 25% FeS, and 12%FeO. 
Such kind of rocks' mass fraction of Si, Mg, and 
Fe is, respectively, about 0.5, 0.62, and 1-05 time s 
that of the bulk Earth (jMcDonough fc Sunl I1995T ). 
For H/He envelopes we use the interpolated h y droge n 
and helium EOS developed by iSaumon et al.1 (|1995| ). 
For wat er we use H2O-REOS, which was applied to 
Jupiter (iNettelmann et all [2008). Uranus and Neptune 
([Fortnev et al.ll2010bl) and in a sligh tly modified version 
to CoRoT-7b (|Valencia et aD 12010 ) . This water EOS 
comprises various water EOS appropriate for different 
pressure-temperature r egimes. It includes th e melting 
curve and phase Ice I (iFeistel fc Wager! I2006D, the sat - 
uration curve and liquid water ([Wagner fc Prufj 120021) . 
vapor and supercritica l molecular water (SESAME 7150, 
iLvon fc Johnsor3[T992]) . and for T > 1000 K and p > 2 g 
cm~ 3 supercritical molecular water, ionic water, supcri- 
onic water, plas ma, ice XII, and ice X based on FT-DFT- 
MD simulations ([French et al.ll2009[ ) . At pressures below 
0.1 GPa, and/or temperatur es below 1000 K, H 2 Q-R EOS 
relies on Sesame EOS 7150 ([Lvon fc Johnsonl[l992[ ). For 
mixtures of hy drogen, helium, and wate r we use H-REOS 
for hydrogen (|Nettelmann et al.l 120081) H9O -REOS for 
water, and an improved version (lKerlevj|2004D o f the he- 
lium Sesame EOS 5761 ([Lvon fc Johnsonlll992t ). Other 
materials are not considered here0. 

4 The label metals comprises all elements heavier than H and 
He. 

5 We do not consider the lighter and more volatile ices CH4 and 
NH3. We performed simple tests by perturbing warm H2O-REOS 
adiabat s with the zero-temperature P — p relations for CH4 and 
NH3 bv lZharkov fc Trubitsvnl 1)19781 ). In a solar C:N:0 proportion, 
this reduces the density by ~ 0.5 g cm -3 if P > 20 GPa, but 
enhances the density for smaller pressures, due to neglcction of 
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2.3. Mass loss 

Mass escape caused by stellar energy input is known 
to occur from the highl y irradiated atmosphere of the 
hot Jupiter HP 209458b (iVidal-Madiar et al l[2003] lYeTIel 



[2001 lErkaev et al.l l2007t iMurrav-Clav et al.l 120091 ) . and 

shown to have a significant impact on the current com- 
position of sup er-Earth CoRoT-7b ([Valencia et al.ll20Tot 
I Jackson et all 120101 ) . While the sun-like star CoRoT- 
7 irradiates the planet with a present X-ray and ultra- 
yiolet energy flux (X UV) F X uv = 5 x 10 5 erg cm -2 s" 1 
(|Valencia et al.ll2010D. GJ 1214 is suppos ed to be an in- 
active M star ([Charbonneau et al.l 120091 ) . Assuming it 
obeys the empirical relation for the surface energy flux 
of M stars, F* xuv /F* boi = (1 - 100) X 10~ 5 (jScalo et al.l 
l200l . then with = 0.00332,0, i?* = O.211i? , the en- 
ergy flux Fxuv = ^ iX uv(-Rp/4a 2 ) received by GJ 1214b 



is only (0.8—80) erg cm 



4) 

i.e. 



(0.16 - 16) x 10~ 



that of CoRoT-7b, at the current time. We then expect 
XUV irradation to have a comparatively lesser, but still 
important, influence on the atmospheric mass through 
time. 

We use t he e nergy limited escape model of 
lErkaev et al.l (j2007f ) to investigate the mass-loss history 
of the planet. With a heat absorption efficiency e = 0.4, 
and a correction factor -Ktide = 0.95 accounting for the 
height decrease of the Roch e-lobe boundary through tidal 
effects (|Erkaev et al.ll2007D . the energy-limited mass es- 
cape rate M esc = (e Fxuv-Rp) (G M p K^)' 1 (following 
I Valencia et aT1l20ToT ) is (0.012 - 2.8) x 10 s g s _1 . This 
is w (0.02 — 6) x 10 -3 the mass loss of present (rocky) 
CoRoT-7b, with a value of (1 - 100) x 1.84 x 10 6 g s" 1 , 
for the fiducial GJ 1214b values M p = 6.55 M ffi and 
R n — 2.678-R®. Correcting M esc ~ r\R p for the altitude 
> R p where the XUV flux is absorbed (jLammer et al.l 



120031 ). the actual value can further rise by a factor of 
10. Within 1 Gyr, this mass loss accummulates to 
(1 - 100) x 10~ 4 M @ . For a low -mass atmosphere of 
only about 1% of GJ 1214b's total mass (which is quite 
possible for a small mean molecular weight atmosphere: 
see § 13. lj) . the fraction of the atmosphere lost during the 
3 to 10 Gyr lifetime of this planet can be large enough to 
influence its cooling behavior. In particular, according 
to the estimates above and assuming constant mass loss 
over time, the initial atmosphere can have been larger 
by (1 — 100) x (0.5 — 1.5)%. We therefore must include 
mass loss in our evolution calculations of H/Hc envelope 
models, as described below. 

2.4. Interior structure and evolution modeling 

Three classes. We consider three classes of hydrostatic 
two-layer interior models of present GJ 1214b assum- 
ing a homogeneous envelope above a rocky core. The 
three classes differ in the materials constituting the en- 
velope. Class I models have a H/He envelope with 
solar He mass fraction Y = Mn e /(Mn + Mhc) 
0.27. Class II models have a pure water envelope ("wa- 
ter worlds"), and class III models a H/Hc/H^O enve- 
lope with variable water mass fraction Z\, and also 
Y = 0.27. The rock core mass of present GJ 1214b 

finite temperature effects which are important for supercritical ices. 
Since these EOS are not of comparable quality to those used for 
all other components, we do not include them here. 



models is found by the condition to match the radius 
R p = 2.678i?© (recently confi rmed through bo th opti- 
cal and infrared photometry dSada et al.l 12010ft ) for a 
given planet mass M p = 6.55M ffi and surface thermal 
boundary condition. We do not consider different possi- 
ble Rp — Mp pairs within the observational error bars 
o~m p = ±1.0M ffi and <tr p = ±0.13i?© as such work 
has already been presented by iRogers k Seageil (|2010|) . 
For the thermal structure of the atmosphere we ap- 
ply the solar composition model atmosphere between 
20 mbar and 10 bar to model classes I and III, and 
to model class II the wa ter model atmosphere betwee n 
20 mbar and 1 bar from iMiller-Ricci k Fortnevl (|2010ft . 
At higher pressure in the radiative atmosphere, we as- 
sume an isothermal temperature — - a reasonable assump- 
tion (Miller-R icci k Fortnevl 120101 ) and given and our 
general understanding of highly irradiated atmospheres. 

Calculating the structure. We choose 20 mbar as the low- 
pressure boundary of our models. This is a choice of 
convenience, since our water EOS ends at this pressure, 
but it is also realistic. The wide-band optical transit ra- 
dius for the planet is at ~10 mbar (E. Miller-Ricci, per- 
sonal communication). This is consistent with the clo ud- 
free atmosphere calculations of iFortnev et al.l (|2003l ) for 
HD 209458b, as well. 

The high-pressure boundaries of the model atmo- 
spheres of respectively 1 and 10 bar are chosen within 
the isothermal part of the atmosphere, before it tran- 
sitions to the adiabatic interior at some pressure P a d- 
For present time (to) structure models, we consider the 
transition pressure P a d(to) & variable parameter and in- 
vestigate the response of the core mass and the cooling 
time on the choice of P a d(to)- 

Since the model atmospheres predict almost constant 
equilibrium abundances of H, He, and H2O, we derive the 
mass density p(P,T(P)) in the atmosphere T(P) from 
an EOS table for constant composition. Given M p , R p , 
and the P — p relations according to the EOS p(P, T(P)) 
in the atmosphere, and constant s(P,T, p(P,T)) in the 
adiabatic interior, we obtain internal profiles m(r), T(r), 
and P(r) by integrating the equation of hydrostatic equi- 
librium dP/dr — —Gmp/r 2 from the surface toward the 
center. Mass conservation M„ 



J Q Rp dr 4nr 2 p(r) 



is en- 



sured by the proper choice of the rock core mass M core . 

The rock core assumption. We assume the rocky core to 
be appropriatly described by an EOS of homogeneous, 
adiabatic 'rocks' at all times (compare Rogers and Sea- 
ger 2010: differentiation into an iron core and a silicate 
layer of Mgr^iFer^iSiOa at uniform temperature). Dif- 
ferentiation will have occurred when temperatures in the 
primitive rock core rose above the melting temperature 
of iron during formation, and will have affected the ther- 
mal evolution. In the Earth, solidification of the inner 
iron core still causes a bouyancy of light elements driving 
convection of the outer iron core, and supports -together 
with gravitational energy release from core shrinking- 
subsolidus convection of the silicate mantle. If the melt- 
ing line of iron rises steepl y with temperature a s indi- 
cated by ab initio data (see iValencia et al"1l2010l for an 
overview), also the central part of a several earth mass 
core of GJ 1214b might transition from liquid to solid 
iron due to high pressure up to ~ 15 Mbar, supporting 
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the assumption of an adiabatic interior. 

Due to the poorly constrained iron mass fraction of 
GJ 1214b and uncertainties in the iron melting line, the 
deep interior could potentially be fully liquid or solid and 
isothermal. Since the equations of state of rocky materi- 
als at high pressure above few Mbar are not well known 
and the e ffect of temperatur e on the P — p relation is 
negligible (jSeager et al.ll2007l ). we believe the P — p rela- 
tion of the rock-EOS used is appropriate for our purpose 
of determining the core mass and its contribution to the 
cooling time. We denote by T core the temperature at 
the core-mantle boundary and assume its time deriva- 
tive T core /dt to be representative for the whole core. 

Calculating the evolution. Chosing P a d(to) for class I and 
II models, and P a d{to) an d Z\ for class III models 
uniquely defines the core mass of resulting interior mod- 
els. A selection of six such present time interior models is 
shown in Table[TJ For class I and II models, we calculate 
the cooling curve by first generating ~ 50 profiles with 
decreasing transition pressures P a d(to) > Pad > 20 mbar 
values, thereby increasingly warmer interiors. For each 
of these intermediate profiles, core mass and composi- 
tion are conserved. As the interior becomes warmer with 
decreasing P a d, the planet radius R p rises. In order to 
obtain the cooling curve R p (t), we integrate the energy 
balance equation 



dE-„ 



dt 



(1) 



backward in time, starting with the present time struc- 
ture models. In Eq.(JTJ, L e ff = 47r R^aT^ s is the net 
luminosity the planet radiates into space, and L cq = 
4nRpaTg q is the stellar energy absorbed. The difference 
L e ff — L cq <~ Ttt is provided by our atmosphere grid 
^int(ff)S), see § 12.1.1 and sets the intrinsic luminosity the 
planet can radiate away from the interior through its at- 
mosphere at a given gravity and internal entropy. Given 
Tint, we can then derive the time interval dt it needs to 
lose the intrinsic energy dE ultl 



dE in 
dt 



dm 



Tds 
~dT 



C, M r , 



dT r , 



dt 



(2) 



Expression ([2]) accounts for the heat loss Sq = T(m)ds 
of each envelope mass shell dm, the heat loss of the core 
due to cooling, and the energy gain L ra dio of rocky core 
material due to decay of radioactive elements (see below). 

Experimental data for the specific heat c v of warm, 
compressed rocks at P > 2 Mbar arc not available. 
Ab initio calculations for iron at a few Mbar and sev- 
eral thousand K (Earth's core conditions) predict c v ps 
0.5 JK-ig" 1 (lAlfe et al. || 2001h. while c„ = 1 .0 J K^g" 1 
was formerly appli ed by Guillot et al.l (119951 ) to the core 



of Jupiter and by iValencia et al.l ( 20101) to the silicate- 
iron interior of CoRoT-7b. We aim to bracket the uncer- 
tainty in c v by chosing c v = 0.5 — 1.0 J Kg. For given 
Ti n t, the time interval required to lose Ei nt rises with c v . 



genie heat. As we will see, our class I and III planet 
models have large rocky cores. Modeling the radiogenic 
heat from the rocky portion of the interior is impor- 
tant to accurately calculate the thermal evolution. For 
Lradio we consider the isotopes 238 U, 235 U, 232 Th, and 



40 K. We adopt isotopic abundances and element abun- 
dances of meteo rites for the elements U, Th, K, and 
Sfl as given in lAnders fc Grevessd (|1989f ). With re- 
spective half-life times in Gyrs of 4.468, 0.704, 14.05, 
and 1.27, and respective decay energies in MeV of 4.27, 
4.679, 4.083, and 1.320, we find a radioactive energy re- 
lease for 1M® of meteoric material on Earth today of 
2.3 x 10 13 Js-yiMe (=: L radi0tmeteor {to)) and 4.56 Gyr 
ago of 2.7 x 10 14 Js-yiMe (=• Lr adi( , >meteor (0)). The 
dominant contribution during this time interval is mostly 
due to (3~ decay of 40 K into 40 Ca. Extension of the ra- 
dioactive decay law to 10 Gyr ago would increase L ra dio 
to 20 x L radio, meteor (0), an unrcalistically high value for 
rock material in the young universe. Of course, for small 
core masses (0.1M p ), cooling of the core does not signifi- 
cantly affect the model's cooling time, but for large core 
mass models it does, and in particular the choice of the 
initial value L ra dio{0) matters. Therefore, we define a 
cosmological luminosity L cosmo := L ra dio,meteor(0) that 
we require each cooling model of GJ 1214b to start with, 
within 10%. Models with long cooling times (10 Gyr) will 
then have L ra dio(to) ~ O.lx L ra dio,meteor(to), and models 
with short cooling times (3 Gyr) will have L ra dio{to) ~ 
2 x Lradto,meteor{to). By this choice of L cosmo wc avoid 
extremely high or low initial values. 

Evolution with mass loss. For a given planet radius, low 
mean molecular weight atmospheres or envelopes are 
also of low mass, and hence such envelopes may lose 
a larger relative mass fraction than large mean molec- 
ular weight atmospheres. We take into account mass 
loss only for the H/He envelope models, e.g. models 
Ia,b in Table [TJ For each neighbored pair of interior 
profiles as given by P a d, we calculate the mass dM 
lost during a time interval dt self-consistently using a 
Newton-Raphson scheme to find the root of the function 
f(dt(dM)) = dM - M esc dt(dM) = 0. 

Other contributions. We did not include the effect of tidal 
heating due to tides raised on the planet by the star. 
This additional amount of energy would tend to delay the 
cooling time, in particular in the past when the eccentric- 
ity e, and hence the tidal heating dE/dt ~ (e 2 /ap) (to 
second order in e, e.g.; iBatvgin et aLll2009b[ ). was larger 
than today. On the other hand, inclusion of tidal mi- 
gration d p of the planet would have reduced the amount 
of ste llar irradition received at early times (jMiller et al.l 
I2009D . A proper treatment of both effects is beyond 
the scope of this pa per since the orbital ecce ntricity is 
not well-constrained (|Charbonneau et al.ll2009() . Instead, 
we make the standard assumption that the formation of 
GJ 1214b yielded an initial heat of accretion that was not 
released before the planet arrived at its current location. 

The cooling time of our thermal evolution models is the 
time between the present state and the time when the 
derivative dT c g/dt approaches — oo, which corresponds 
to a hot start that is insensitive to the initial condi- 
tions. Whi le such a treatment of young planets is com- 
mon (e.g.; IBaraffe et aLll2003[) . it actually ignores the 

6 The Si mass abundance of meteorites (10.65%) is unequal that 
of the rock EOS (16.8%). 

7 These values are in part an averag e over the surpri singly wide 
spread of literature values, see e.g. http://ic.lbl.gov/education/ 
and in part over decay chains ( 4t) K) 
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Table 1 

GJ 1214b structure model input and output data 



Label 




Pad{t ) 




envelope material 


M 


k 2 


Tint 


age 




(AT©) 


(bar) 






(10 7 g/s) 




(K) 


(Gyr) 


la 


6.464 


300 





H/He 


1.84 


0.0183 


42.7 


3.1 


lb 


6.434 


800 





H/He 


1.84 


0.0170 
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process of planet formation which has been shown to 
alter t he luminosity durin g the first tens of millions of 
years (|Fortnev et al.l [20051 ) . This is important for mass 
determinations from cooling tracks, whereas in our case 
it just might induce an error of ~ 10 Myrs to the calcu- 
lated cooling time. 

Love number hi. We follow the call bv lRagozzine fc WolJ 
(|2009f ) to tabulate the tidal Love Number ki values of 
representative models. This quantity is a planetary prop- 
erty which solely depends on the internal density distri- 
bution. If known, ki imposes an additional constraint 
on interior structure models. Physically, ki quantifies the 
quadrupolic gravity field deformation at the surface in re- 
sponse to an external perturbing body of mass M, which 
can be the parent star, another planet, or a satellite. 
M causes a tide-raising potential W(r) = J2 n =2 W n = 
(GM/a) ^ n=2 (r/a)™P n (cos?/>), where a is the distance 
between the centers of mass (planet and body), r is the 
radial coordinate of the point under consideration in- 
side the planet, i\) is the angle between the planetary 
mass element at r and the center of mass of M at a, 
and P n are Legendre-Polynomials. Each external po- 
tential's pole moment W n {r) induces a change V n (r) in 
the corresponding degree of the planet's potential, where 
V n (Rp) = k n W n (R P ) defines the Love numbers. For 
th e calculation of ki we follow the approach described 
in lZharkov fc Trubitsvnl (fi97l) . 

To first order in the expansion of the planet's potential, 
ki i s proportional to the gravitational moment Ji (see 
e.g. lHubbardlllQM §4). That is why measuring ki will 
provide a constraint for extrasolar planets that is equiv- 
alent to J2 for the solar system planets. In particular, ki 
is known to be a measure of the central condensation of 
an object, which can be parameterized by the core mass 
within a two-layer model approach. However, as for Ji, 
the inverse problem — the deduction of the internal den- 
sity distribution of the planet from ki, is non-unique. 

We stress that the Love number ki is a 'poten- 
tially observable parameter if the planet's eccentric- 
ity is non-zero. It can be obtained with the help of 
transit light curves that contain information about the 
tidally induced apsidal precession of close-in planets 
(|Ragozzine fc W olf 2009), or for specific two-planet sys- 
tems can al so be derived from measuring the orbital 
parameters (|Batvgin et al.ll20"09b( ). given that the plan- 
ets are in apsi dal alignment and on co-planar orbits 
(IMardlindl2010h . 

3. RESULTS 

The framework of the three structure classes I-III is 
used to investigate the possible set of models with respect 



to core mass, internal pressures and temperatures, water 
to rock ratio, and the H/He mass fraction. 

3.1. GJ 1214b with a H/He atmosphere and a rock core 

For our H/He envelope+rock core models we find a nar- 
row core mass fraction range of 0.975—0.995 for a wide 
pressure range 50 < P ad (to) < 4000 bar (Fig.QJi). De- 
spite its low mass fraction, the H/He envelope extends 
over 0.37 R p (Fig.[Tb) independent of the envelope mass. 
This is because of its high temperature, which increases 
from 1030 K at P ad {t ) to 2400-4200 K at - 35-90 kbar 
at the envelope-core boundary. At these conditions, hy- 
drogen is molecular throughout the envelope according 
to the SCvH-i EOS. 

Not all of the models shown in Fig. QJ,,b are consis- 
tent with an age t+ = 3 to 10 Gyr of the star GJ 1214. If 
P a d(to) decreases (increases), internal entropy rises (falls) 
and the planet will need less (more) time to cool down 
to this state. This behavior is illustrated by the cooling 
curves in Fig. [Tp, according to which an age of 3 Gyr re- 
quires P a d(to) > 300 bar; and a much longer cooling time 
of 7.2 Gyr is obtained for P a d(to) = 800 bar, whereas 
Pad(to) = 100 bar would give a cooling time below 1 Gyr. 
With a mass loss rate according to Fxuv / Fbo\ = 10 -5 as 
assumed for the cooling curves in Fig.[T};, a cooling time 
of 10 Gyr can not be obtained through a further increase 
of P a d beyond 800 bar if L ra( no(0) is not to drop be- 
low 0.9 x L cosrno . The colder the interior, the lower T- mt 
as predicted by our model atmsphere grid, and hence 
the intrinsic energy then can be transported through the 
radiative atmosphere. Enhancing P a d(to) from 300 to 
800 bar lowers T int from 42.7 to 31.8 K (Tabled}. For 
even colder interiors, the atmosphere is no longer capa- 
ble of radiating away the heat generated by radioactive 
decay, which would contradict the assumption of such a 
cold interior. 

The cooling time increases with the mass loss rate 
(FigEDi)- For stellar XUV radiati on F X uv /F bol = K T 5 - 
10~ 4 as typical for quiet M-stars (|Scalo et al.ll2007D . this 
enhancement is small (see Fig. [T}i) and we obtain es- 
sentially the same range of P a d(to) arL d hence structure 
models. Of GJ 1214b's initial total mass (initial H/He 
envelope mass), only 0.0054-0.013% (0.41-0.74%) is 
lost if Fxuv/^boi = 10" 5 , and about 0.1% (6%) if 
Fxvv/Fboi — 10~ 4 . On the other hand, for a perma- 
nent, strong irradiation i*xuv/-Pboi = 10 -3 as observed 
for young, active M-Stars listed in the ROSAT catalogue, 
GJ 1214b would have lost 1.8% (53%). These numbers 
are in agreement with the rough estimates in §12.31 In the 
last case, existence of a thin H/He atmosphere to-date 
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Figure 1. Cooling curves and internal structure of models with H/He envclope+rock core. All models are adjusted to have R p (to) = 
2.678-R® by the choice of the core mass, and M p (to) = 6.55M0 at present time to- Upper panels: Mass (panel a) and radius (panel b) of 
the core (solid line) and where the atmosphere transitions into the adiabatic interior (dashed) of present time models with 50 < P a d(to) < 
4000 bar. Numbers are pressure in kbar and temperature in K at the envelope-core boundary for selected models. Single interior models run 
vertically. The gray shaded areas are a guide to the eye for the allowed P a d(to) range that is consistent with the thermal evolution. Lower 
panels: Evolution of radius with time. Panel c: its dependence on P ac ;(io)> panel d: on the insolation as parameterized by i^xuv/^bob 
dashed: model la, dotted: model lb (see Table[T]l, dot-dot-dashed: an intermediate one with P a d(to) = 500 bar. These cooling curves are 
for a specific heat of the core material c v = 1.0 J/gK and an insolation i^xuv/^bol = 10 -5 (black), 10 — 4 (thin solid gray), and 10 — 3 (thick 
solid gray). Circles in panel c indicate the profiles during evolution when P a d(t) = 100 bar. The boxes indicate the observational error 
bars of radius (vertical extension) and stellar age (horizontal extension). 



becomes less likely since it begins to require fine-tuning 
of the initially accreted H/Hc envelope mass. We have 
found cooling tracks with an age of 10 Gyrs or more 
only if P > 800 bar and the mass loss rate is high, or 
L ra dio(t) <C L cosmo , which we do not favor. 

Models with 300 < P a d < 800 bar as constrained by 
our evolution calculations have M core = 6.434 — 6.464M® 
and k 2 = 0.0170 — 0.0183 implying a high degree of cen- 
tral condensation. The presence of an iron core could 
even enhance the central condensation. Therefore, al- 
though the H/He layer is low in mass, it significantly 
strengthens the property of central condensation com- 
pared to a closer to zero-mass atmosphere planetary ob- 
ject such as the Earth , the theoretical k 2 value of which 
is - 0.3 (|Zhandll99ll ). 

Class I models are closest to giant planets that formed 
within the snowline of the disk and did not have enough 
time and/or material in their surrounding to accrete a 
massive H/He envelope. Class I stru cture models can 
best be compared to 'case F models by iRogers fc Seagerl 



(|2010f ) . where the difference in composition assumptions 
can be reduced to the core (undifferentiated rock core 
versus differentiated iron-silicatc-water-ice core in their 
models), and a slightly different envelope He abundance 
of respectively 0.27 and 0.28. 

Our obtained H/He mass fraction range, 1.3—1.8%, is 
due to the uncertainty in T; nt . This range is much smaller 
than theirs (9 x 10 -3 — 6.8%), which includes uncertain- 
ties from the la errors of M p and R p contributing an 
uncertainty of 0.4—1.6%, from their Tj nt values used con- 
tributing up to 1.5%, and from the uncertainty in core 
composition (pure iron, iron-silicates-water, or pure wa- 
ter). Since iron is included in our rock EOS, we consider 
(1.8%+1.6%=) 3.4% a reliable upper limit of the H/He 
mass fraction if observations reveal a low-mean molecu- 
lar weight atmosphere. Recent observations of GJ1214b's 
atmosphere in the 0.78 to 1 ^tm wavelengths range with 
the VLT facility's UT1 telescope suggest a mean molecu- 
lar weight of 5 g/mol or more and thus di sfavor a H/He- 
dominated atmosphere dBean et al.ll2010 |). On the other 
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hand, the presence of clouds or hazes in such an atmo- 
sphere could mimic a short scaleheight, hence high mean 
moelcular weight, and cannot be excluded by current ob- 
servations and model atmospheres. 

3.2. GJ 1214b as a water planet with a rock core 

The derived core mass of water envelope+rock core 
models responds much more strongly to a change of the 
onset of the adiabatic part of the envelope than models 
with H/Hc envelope. At P a d(to) ~ 400 bar, the core 
mass reaches zero: a pure water planet (Fig. [2^). Deeper 
isothermal regions would only be possible if some amount 
of the pure water planet would be replaced by lighter 
elements such as methane or ammonia. Warming up the 
deep envelope by an outward shifting of the onset of the 
adiabatic region below 50 bar is accompanied by a strong 
rise in core mass. A solar water:rock ratio of 2.5 occurs 
for P ac i = 5 bar — when the deep interior is extremely hot. 
Smaller ratios would be possible only if the isothermal 
region is allowed to disappear. We next consider whether 
such models are consistent with the cooling time. 

Because of the relatively small core mass fraction of 
class II models, varying the specific heat of the core has a 
negligible effect on the cooling time. For 80 < P a d(to) < 
300 bar the resulting cooling times are consistent with 
r*, see Tabled] and Fig.[2j;. Those models have M core = 
0.20 — 0.87 Mq, and a waterirock ratio of 6.5— 31. This is 
much higher than the solar ice to rock ratio (I:R)® ~ 2.5, 
which would give a cooling time shorter than 3 Gyr. 

In contrast to models Ia,b, models IIa,b are weakly cen- 
trally condensed as parametrized by their high k^ values 
of respectively 0.57 and 0.74. 

Clas s II models are different from those by iFu et al.l 
(|2010f ) who consider cold water planets with ice or liquid 
ocean layers above silicate/iron cores. GJ 1214b is not 
that cold. Our models also differ from t he water steam 
atmos phere models of CoRoT-7b by IValencia et al.1 
(|2010f ) who find water envelope + silicate/iron core mod- 
els where water in present CoRoT-7b contributes at most 
10% to the total mass and is in the vapor phase or super- 
critical molecular phase. They also differ from the wa - 
ter steam atmosphere models bv lRogers fc Seagerl (|2010f ) 
who describe the deep interior by a water-ice equation of 
state, whereas according to the phase diagram (Fig. ED ) 
(updated from lFrench et all 120091 : IValencia et al.ll2010¥ . 
water would be in the plasma phase in GJ 1214b. 

In our class II models, water transitions from the vapor 
phase to supercritical molecular water still in the atmo- 
sphere, becomes dissociated into an ionic fluid at about 
4000 K and 0.2 Mbar, and finally fully dissociated and 
ionized for T > 5100 K and P > 1 Mbar (Fig.Eb) form- 
ing a plasma with el ectronic conductivity > 100 fi/cm 2 
(jRedmer et al.|[20TTI ). 

Under these circumstances, a planetary interior can be 
able to maintain a dynamo generating a dipolar magnetic 
field. In Uranus and Neptu ne, the magnetic field may be 
generated in a thin shell ([Stanley fc Bloxhanl 120061) of 
possibly ionic water (|Nellis et al.1119881) . This view of the 
cold (T e g « 59 K) outer planets Uranus and Neptune 
is supported by the Neptune adiabat in Fig. [3] while in 
GJ 1214b, since it is warmer, the fluid conductive en- 
velope would extend down to the small core, somewhat 
akin to Jupiter, and therefore preferably lead to a dipo- 
lar field, according to the field geometry considerations 



bv IStanley fc Bloxhaml (f2006h . 

3.3. GJ 1214b as a H/He/H20 planet with a rock core 

This structure type resembles the outer envelope of 
Uranus and Neptune if Z\ < 0.4, or their inner envelope 
if Z\ > 0.7 ac cording to three-layer Uranus and Nep- 
tune models of iFortnev fc Nettelmannl (|2010f ). A strong 
enrichment in metals of the outer H/He layer in Uranus 
and Neptune is necessary to match the gravity field data, 
and also some admixture of light elements in the deep in- 
terior. 

In § 13.21 we have seen that water envelope+rock core 
models of GJ 1214b yield supcrsolar I:R ratios of more 
than 2.6 x solar. Lowering this ratio can be achieved by 
replacing water with hydrogen and helium. The two lim- 
iting cases of this implementation are a structure where 
a H/He layer is on top of a water layer, or a homoge- 
neous mixture of H/He and H2O. We find that in the 
first case, such differentiated three-layer models (H/Hc, 
water, rock) can not have 1x(I:R)q and be in agree- 
ment with r*, the reason of which is the following. Class 
I and class II models require a radiative atmosphere 
down to 80 bar or more at present in order to meet 
t*. If composed solely of H/He, the atmosphere ex- 
tends over about 0.4 i?®. In order to match a remain- 
ing radius r{M p ) » 2.3 i?®, the core mass fraction of 
the rema ining water+core bo dy is of the order of 20- 
50% fsee IValencialFaII[20Tol . Fig. 8). We find I:R< 0.9 
and M core = 3.5 — 4.2 M® for this case of differentiated 
models. Increasing P a d increases the depth of the thin 
H/He atmosphere, thereby lowering the I:R ratio even 
more. Increasing the planet's mean density within the 
la error bars of M p and R p allows for I:R up to at most 
O.56x(LR) . 

Consequently, the only way to obtain a solar I:R ratio 
is to limit the radius of the H/H e atmosphere by enhanc- 
ing its mean molecular weig ht (jMiller-Ricci et al.|[2009h 
through admixture of water. Here we consider the case 
of equal metallicity in the radiative atmosphere and in 
the adiabatic envelope (our class III models) as parame- 
terized by the water mass fraction Z\ . Figure [Hi shows 
the change of the I:R ratio of single models. The water to 
core mass ratio rises moderately up to Z\ = 0.9, passes 
lx(I:R)® at Z\ = 0.95, and then rises rapidly up to the 
values found for class II models. This behavior depends 
very weakly on the choice of P a d{to)- 

For class III models, the resulting H/Hc mass fraction 
of the planet (see Fig.[3Ji) is about 2—3 times larger than 
in case of a H/Hc layer on top of a water layer. It reaches 
the maximum for Z\ « 0.80 (the core mass must not be 
too large, requiring a high metallicity, and also Z\ not 
too close to 1). We find a planetary H/He mass fraction 
< 5.9% if P ad (to) = 400 bar as in Fig.[§i,b, and slowly 
rising with P a d{to) up to 7% (colder envelopes reduce 
the core mass). However, a colder present time interior 
would take longer than 10 Gyrs to cool. For the cooling 
curve calculations we choose a metallicity Z\ = 0.85, 
which gives a H/He mass fraction close to the maximum 
value but also a core mass below 2/3 M p , hence a real 
alternative to classes I and II. Figure 0J; shows that the 
isothermal region of present GJ 1214b must end between 
120 and 400 bar to give consistency with a cooling time 
of 3 to 10 Gyrs. 

With Z 1 = 0.85, T core = 5730 K, and P core = 
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Figure 2. Cooling curves and internal structure of models with water envelope+rock core. All models are adjusted to have R p (tg) = 
2.678-Rg) by the choice of the core mass. Left (panel a): Core mass (solid line) and mass shell where water in the envelope enters the 
plasma phase (dashed line) of present time structure models for transition pressures 5 < P a d(to) < 400 bar. Single interior models run 
vertically. The gray shaded areas are a guide to the eye for the allowed P a d(to) range that is consistent with the thermal evolution (in 
panel c). Middle (panel b): Same as (a) but radius coordinate. Numbers at models highlighted by filled circles are pressure in Mbar and 
temperature in K. Right (panel c): Evolution of radius of structure models Ha (dashed) and lib (solid), see Table[TJ Circles indicate those 
profiles during evolution when P a d(to) = 50 bar. 
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Figure 3. Water phase diagram for 20 mbar< P < 100 Mbar and 
GJ 1214b water envelope+rock core models with P a d(to) = 100 bar 
(orange dashed- dot) and 300 bar (red d ashed-dot). The Nep tune 
profile (blue dashed-dot) is adopted from [Rcdmcr ct al. (2 0111 V 

1.4 Mbar (Fig. 2b), model Illb resembles the interior of 
Uranus and Neptune in com position and temperature 
(jFortnev fc Nettelmannll2010D . Lower in total mass, the 
pressure does not rise up to 5—7 Mbar as in the outer 
solar system giant planets, so that water will not adopt 
the superionic phase according to the phase diagram of 
water, but remain in a fluid state in GJ 1214b (Fig. [3]). 
This property bolsters our assumption of a homogeneous 
mixture of water with hydrogen and helium. 

4. DISCUSSION 

4.1. Structure assumptions 

Our GJ 1214b interior models rely on a separa- 
tion of the interior into a rock core and one ho- 
mogeneous envelope of the same composition as in 
the visible atmosphere. In contrast, giant and ter- 
restrial planets in the solar system are not success- 
fully described by such a two-layer structure but re- 
quire the assumption of various internal layer bound- 
aries to be consistent with the atmospheric He 
abundance and the gravity field data (giant plan- 



4.1.1. class III: H/He phase separation? 

For the giant planets, a layered structure is suggested 
in part because of a measured atmospheric depletion in 
He indicating H/He phase separation, and in part be- 
cause of relatively low measured J4/ J 2 ratios indicating 
an enhancement with metals in the deep interior. Ac- 
cording to experim ental and theoretical data on H/He 
demixing (see e.g. lLorenzen et al.l 120091 : iMorales et al.l 
I2009T ) . immiscibility of He in H might also occur in 
our class III models close to the core-mantle bound- 
ary. Given the low interior temperatures (relative to 
Saturn) this could very strongly deplete mos t of the 
envelope in helium (|Fortnev fc Hubb ard 2004|). While 
accurate measurements of the atmospheric He abun- 
dances are extraordinary difficult to perform without the 
use of entry probes even in case of solar giant planets, 
this effect can affect the current depth of the isother- 
mal region (120—400 bar) that we have derived from 
the cooling time calculations. This would be a com- 
pletely different pathway towards a helium-depleted hy- 
drogen a tmosphere than the ou tgassing scenario dis- 
cussed in IRogers fc Seagerl (|2010f ) . On the other hand 
if Z\ <C 0.85, then the envelope does not become mas- 
sive and dense enough for He sedimentation to occur. 

4.1.2. class II: Incomplete Differentiation? 

Our class II models require a process that causes a 
downward sedimentation of rocks in order to separate 
out a water layer, as in Ganymede (|Kirk fc Stevensonl 
[19871) . Up to now, there is no experimental evidence 
that water and silicates or iron become immiscible under 
high pressure, and the timescale for gravitational set- 
tling before the onset of convection, terminating grav- 
itational settling, is essentially unknown for Ganymede 
(jKimura et al.l [2009D . Since the larger primordial heat 
deposited in the ~ 260 times more massive GJ 1214b 
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Figure 4. Cooling curves and internal structure of H/HC/H2O cnvelopc+rock core models with various envelope metallicities. Left: The 
dependence of the core mass (in M p ) on the envelope metallicity Z\ of present time structure models with P a d(to) = 400 bar. Numbers 
at selected models are the I:R ratio (first number) and the entire planet's H/He mass fraction in percent. Model Illb is highlighted (filled 
circle). Middle: Core radius in R p . Numbers are pressure in Mbar and temperature in K at the envelope-core boundary for the same 
models as in the left panel. Flight: Evolution of radius with time for models Ilia (dashed) and Mb (solid), see Table[JJ using specific heat 
of the core material c v = 1.0, or 0.5 J/gK (dotted). Circles in panel c indicate those profiles during evolution when P a d(t) = 100 bar. 



might have caused a rapid onset of convection, a wa- 
tcr+rock interior of GJ 1214b is possibly not fully dif- 
ferentiated. Internal layer boundaries dividing regions 
with different water to rock ratios can not be excluded, 
and the core rather be an ice- rock mixture as suggested 
for Callisto (INagel c t al. 20QH) than pure rocks. On the 
other hand, more detailed envelope models of GJ 1214b 
would be undcrdetermined by current observational pa- 
rameters. 

4.1.3. class I: Choked off giant planet formation? 

The core mass of our class I models is within the 
range of rock core masses cur rently proposed for Jupiter 
(|Fortnev fe Nettelmannll2010[ ). With an initial H/He at- 
mosphere of only 1—2% M p , GJ 1214b appears to be a 
giant planet whose envelope mass accretion was choked 
off during formation. It is of general interest for more 
super-Earth planets to be detected to see whether this is 
a common outcome of planet formation. 

4.2. Composition 

From interior and atmosphere models that arc con- 
sistent with the observationally derived parameters T Gq , 
M p , and R p after 3—10 Gyrs of cooling, the most cer- 
tain conclusion we can draw about the composition of 
GJ 1214b is a metallicity of 94-100%. The lower limit 
can further shrink somewhat if the planet in reality is 
dry, and the mass fraction of water in our class III models 
then resembles a mixture of H/He and rocks. In contrast, 
the mass fractions of water, used as a proxy for the ices 
H2O, CH4, NH3, and H2S, is essentially unconstrained 
(0—97%) as is the mass fraction of rocks (3—99%). 

Two further observables we can hope to attain in the 
near future are the Love number k% from transit timing 
variations or from the sh ape of the transit light curves 
(Ragozzinc & Wolf 2009), and second the mean molec- 
ular weight in the atmosphere from transmission spec- 
troscopy, in particular from the wavelength dependence 
of star light absorption in the planetary a tmosphere dur- 
ing transit (jMiller-Ricci fe Fortnevl[20ia) . 

4.2.1. Love number k-y 
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Figure 5. Collection of {M CO re, &2} pairs obtained for class I- 
III models of GJ 1214b (diamonds), that meet the age constraint 
(filled) or not (open); and measured &2 values with error bars 
of Moon, Mars, and Earth (M M E) (circles). According to 
the definition of M core used for GJ 1214b, M M E would have 
Mcore > 0.999999M P , which we placed at 0.99Af p . Considering al- 
ternatively M core as the iron core of these objects shifts the filled 
circles to the open circles. This figure demonstrates degeneracy of 
ki in two-layer models with thin atmospheres. 

The &2 values of models from our classes I-III differ 
greatly from each other and thus we consider ki a useful 
quantity to discriminate between atmospheres of differ- 
ent mean molecular weight. The trend of decreasing k% 
value with increasing core mass that is known for two- 
layer models of J upiter-mass giant pla nets with a core 
and one envelope (jBatvgin et al.f 2009b). o r for n =l poly- 
tropic planets in general (jKramm et al.1 120 111 ) , is con- 
firmed by our two-layer models of GJ 1214b. An illustra- 
tion of the k2(M core ) behavior is shown in Fig.0 which 
contains the same data as TableTJJand a l so some interme- 
diate points as well as two solutions of the H/He envelope 
models with respectively P a d{U)) = 1 and 2 kbar, which 
are too cold as explained in § 13.11 

However, when the H/He atmosphere becomes thin 
and its mass low enough (< 2%), this planet begins 
to more closely resemble a relatively homogeneous rock 
body than a core+envelope planet, and ki rises again. 
This degeneracy of k^ that is otherwise well-known for 
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multi-layer models is -at first glance- a surprising find- 
ing in two-layer super-Earth models. As Fig.[5] suggests, 
pure H atmosphere models of GJ1214b (|Rogers fc Seagerl 
|2010| ) would approach the k^ values of solar system bod- 
ies such as the Moon, i.e. higher than that of H/He at- 
mospheres. 

Moon, Mars, and Earth can be considered as rem- 
nant protoplanetary cores (M core = M p ) that did not 
accumulate sufficient mass to accrete a significant enve- 
lope. On th e other hand, be cause measured k^ values 
of Earth (seelRav et al.l l2001IL Mars (lYoder et all 120031: 
iKonopliv et al.ll2006f l. and Moon ((Zhang|ll992j) cannot be 
explained by a homogeneous rock interior but require the 
assum ption of a dense, iron-rich core (e.g. lZharkov et all 
f2009h . they can also be considered as objects with small 
(iron) core and large (silicate) mantle (M core <C M p ). 
This places Mars and Earth close to our GJ 1214b mod- 
els with H/Hc/Z envelope in the ki — M core diagram, 
implying that ki is a degenerate quantity with respect 
to composition, too. Separation into layers of differ- 
ent composition owi ng to phase differentia tion and phase 
transitions in rocks ((Valencia et alJl2007| ) has likely also 
occured in GJ 1214b. Such an advanced treatment of 
a core would enhance the level of central condensation, 
and result in even lower minimum ki values below 0.017. 
The same trend is expected for inclusion of solid-body 
effects, which will be important for colder super-Earth 
planets and those with less extended, less massive at- 
mosphere. Theoretical ki values for the terrestrial ob- 
jects are in good agreement with the observations when 
making the simplifying assump tion of an elastic interior 
(lZhanelll992l: lYoder et al.l[2003h (in which case the shear 
modulus becomes frequency-independent, i.e. ki a static 
Love number), while their observed moments of inertia 
are close to 0.4 indicating a nearly homogeneous interior, 
the theoretical Love number ki of which would approach 
1.5 if the body were fluid and compressible. 

4.2.2. Mean molecular weight 

Class III models with lx solar I:R ratio are not ex- 
cluded by our thermal evolution calculations. Those 
models require Z\ = 0.95, corresponding to a mean 
molecular weight [i = 13.4 gmol^ 1 and 440xsolar 0:H 
particle number ratio (assuming a solar 0:H of 0.85 x 
10~ 3 according to I Anders fc Grevessd 119891 ). However, 
we do not consider such models a realistic description 
of GJ 1214b. For rapid translations between metallicity 
and atmospheric /x values from transmission spectroscopy 
observations, we present the simple relations between /x, 
0:H, and Z x in Fig.© 

In particular, our class III models with Z\ = 0.85 have 
[i = 8.9gmol _1 , a 270x solar atmospheric 0:H ratio, 0.2 
times solar bulk I:R ratio, and ~ 5.8% of M p is H/Hc. 
For Z\ < 0.55, our class III models become severely rock- 
dominated with M core > 0.9Aip, /1 < 4.4gmol _1 and 
0:H < 92 x solar. A 50 x solar 0:H ratio as of our model 
atmosphere grid applied to the thermal evolution calcu- 
lations (see § I2.1|) implies Z\ = 0.38. For /1 = 5gmol _1 
as suggest ed from the first tr ansmission spectrum obser- 
vations bv iBean et all (|2010T ). Fig. [6] gives Z\ = 0.7 and 
according to Fig.|H GJ 1214b could be mostly rocky with 
an I:R ratio of 0.15, where the H/HC/H2O atmosphere 
contributes ~ 18% to the planet's mass and 40% to its 
radius. 




i i 
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Figure 6. Mean molecular weight {dashed) and 0:H ratio (solid) 
of H/HC/H2O mixtures with water mass fraction Z and Y=0.27. 

5. CONCLUSIONS 

Our results for the composition of the super-Earth 
mass planet GJ 1214b confirm that it has a gaseous at- 
mosphere atop a fluid envelope. We find a minimal total 
H/Hc mass fraction of 1.3% for pure H/He envelopes, 
which can rise to 5—6% if the envelope contains 60—90% 
H2O in mass, and even further if also silicates and iron 
are mixed into the envelope. 

Water in GJ 1214b does not solidify within 10 Gyrs of 
cooling, and it not a liquid as is found on Earth's surface, 
but becomes a plasma if its abundance is high {Z\ > 0.8). 
This leads to a large water mass (> l/3M p ), where the 
deep internal matter is warm (T > 5500 K) and dense 
(P > 1 Mbar). 

The intrinsic heat loss of GJ 1214b after 3—10 Gyr of 
cooling corresponds to Tj nt = 32 — 52 K, where the lower 
bound depends on assumptions about the heat produc- 
tion by radioactive elements and slightly on mass loss. 
By our self-consistent mass-loss calculations we conclude 
that GJ 1214b is a genuine super-Earth that has lost an 
insignificant amount of its initial mass, unlike CoRoT-7b. 

Nevertheless, some of our computed models may fall 
into the realm of the unlikely. These include class I mod- 
els that have lost more than 50% of their initial H/Hc en- 
velope (i.e. 2%M p lost) due to extraordinary XUV irra- 
diation and thus require fine-tuning of their initial H/Hc 
mass fraction. Moreover, extremely low-mass H/He at- 
mospheres in general might be an unlikely outcome of 
planet formation. The unlikelihood of forming a massive 
planet with an extreme ice-to-rock ratio of more than 6 
casts our class II ("water world") models in doubt. We 
instead favor what is thought to be the most Uranus- 
and Neptune-like planet models, class III. It includes en- 
velopes with ~ 85% water by mass mixed into H and He 
atop a rock core with about 0.2 x solar bulk icerrock ra- 
tio. These models have higher H/He masses (5.8% of the 
planet mass) than class I models with pure H/He atmo- 
spheres, lessening worries about the survival of a H/Hc 
envelope in the face of mass loss. A general outcome 
of these models is T int — 42 ± 6 K, while the condi- 
tions at the core-mantle boundary arc P ~ 1 Mbar and 
T w 5700 K. The atmosphere can be strongly depleted in 
helium due to H/He phase separation deeper inside the 
planet (if Z\ ~ 0.85) or not be depleted in He if Z\ is 
much smaller than this. 

Calculated Love numbers of water envelope models 
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that are consistent with the observables M p , R p , T c<i , 
and the age, are in the range ki = 0.58—0.74, whereas 
pure H/He envelope models have &2 ~ 0.018, and our fa- 
vorite models ki ~ 0.15. An observational determination 
of &2 and the atmospheric mean molecular weight is cru- 
cial for determining the envelope metallicity and the core 
mass of this and other planets that are located along the 
M — R relation of water planets. However, ambiguities 
in composition are an inherent property of such planets. 

We thank E. Kcmpton (formerly Miller-Ricci) for dis- 
cussion and the delivery of her solar composition and 
water model atmospheres, T. Guillot and D. Valencia 
for valuable comments on the radioactive luminosity of 
rock material, and F. Sohl for helpful discussions on dy- 
namic Love numbers. This work was supported by the 
NASA, under grants NNX08AU31G and NNX09AC22G, 
and DFG under grants RE 882/11 and RE 882/12. 
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